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Abstract. Results on transverse mass spectra of neutral pions measured at central rapidity are presented for impact 
parameter selected 158-yl GeV Pb-l-Pb, and Pb-l-Nb collisions. The distributions cover the range 0.5 GeV/c^ < 
mo < iGeV/c^. The change of the spectral shape and the multiplicity with centrality is studied in detail. In going 
from p-l-p to semi-peripheral Pb-l-Pb collisions there is a nuclear enhancement increasing with transverse mass similar to 
the well known Cronin effect, while for very central collisions this enhancement appears to be weaker than expected. 
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1 Introduction 



2 Experiment 



Heavy ion reactions at ultrarelativistic energies provide infor- 
mation on nuclear matter at high energy density (for reviews 
see e.g. Hadron production is generally considered to 

be sensitive to the late freeze-out stage of the collision, when 
hadrons decouple from one another. Already from the experi- 
mentally determined shape of the transverse mass spectra it is 
evident that heavy ion reactions are not merely a superposition 
of nucleon-nucleon collisions 

The broadening of the transverse mass spectra in pH-A col- 
lisions compared to pH-p (Cronin-effect [0) has been attributed 
to initial state multiple scattering of partons [^^. However, 
models which attempt to describe nucleus-nucleus reactions, 
like the Monte-Carlo programs VENUS [0 or RQMD 
rely on the assumption that final state rescattering plays an im- 
portant role in determination of the momentum distributions of 
the hadrons. If final state rescattering is very strong, the notion 
of local kinetic equilibrium may be valid, which is the basic as- 
sumption for a hydrodynamical description (see e.g. [ 12 , jj]]). 

It remains unclear, whether equilibration is attained. If it 
occurs, it will be important to isolate the contributions of ther- 
mal, pre-equilibrium, and initial state processes to the hadron 
yield. Since the latter are expected to dominate at large trans- 
verse momenta, it is hoped that systematic studies of hadron 
spectra over a large range in momentum might allow to disen- 
tangle these contributions. In the analysis of central reactions 
of Pb H- Pb at 158 A GeV it is seen that both predictions of per- 
turbative QCD [ p^p^ and hydrodynamical parameterizations 
1 19 1 can describe the measured neutral pion spectra reasonably 
well. It is particularly surprising to observe that, on the one 
hand, a pQCD calculation gives a reasonable description also 
at relatively low momenta, while on the other hand, a hydrody- 
namical parameterization can provide a good description, even 
at very high momenta. 

The understanding of the relative contributions of the var- 
ious soft and hard processes in particle production is espe- 
cially important in view of the recent interest in the energy 
loss of partons in dense matter [|l4 15 1, generally referred to as 
jet quenching, as a possible probe for the quark-gluon-plasma. 
One of the suggested experimental hints of jet quenching is the 
supp ression of particle production at high transverse momenta 
|16|. In order to confirm such an interpretation, it is impor- 
tant to study other possible nuclear modifications of particle 
production in detail. More information in this respect may be 
gathered from the variation of the particle spectra for different 
reaction systems or different centralities. First attempts in this 
respect have been discussed in [[l7|], where the neutral pion av- 
erage pt as a function of the centrality was shown to rise from 
pH-p collisions to peripheral Pb-nPb collisions and to saturate for 
medium central and central Pb-nPb collisions. It was also seen 
that the neutral pion yield for reactions with more than about 30 
participating nucleons Npart exhibits a scaling as N^^^.^ with a 
power a 1.1 which is approximately independent of pr- 

In the present paper we will present a detailed study of neu- 
tral pion transverse mass spectra in the range 0.5GcV/c^ < 
rriT — mo < 4.0GeV/c2 and 2.3 < y < 3.0 for collisions of 
PbH-Pb and Pb-nNb at 158-^ GeV for different centraUties. 



The CERN experiment WA98 [ J20| ] is a general-purpose appara- 
tus which consists of large acceptance photon and hadron spec- 
trometers together with several other large acceptance devices 
which allow to measure various global variables on an event- 
by-event basis. The experiment took data with the 158-v4 GeV 
208pb beams from the SPS in 1994, 1995, and 1996. The data 
presented here were taken during the 1995 and 1996 beamtimes 
at the CERN SPS. The layout of the WA98 experiment as it ex- 
isted during the final WA98 run period in 1996 is shown in 
Fig.0. 

Neutral pions are reconstructed on a statistical basis from 
their two-photon decay, using photons measured with the 
LEDA spectrometer in the pseudorapidity interval 2.3 < 77 < 
3.0. This detector is located 21.5 m from the target and consists 
of 10,080 modules. Each module is a4x4x40 cm'^ (14.3radia- 
tion lengths) TFl lead-glass block read out by an FEU-84 pho- 
tomultiplier The high voltage is generated on-base with cus- 
tom developed Cockcroft-Walton voltage-multiplier type 
bases which are individually controlled by a VME processor. 
The photomultiplier signals are digitized with a custom-built 
ADC system [p3||. Twenty-four lead-glass modules are epox- 
ied together in an array 6 modules wide by 4 modules high 
to form a super-module. Each super-module has its own cali- 
bration and gain monitoring system based on a set of 3 LEDs 
and a PIN-photodiode mounted inside a sealed reflecting front 
cover dome p4| | . The photon spectrometer is separated into two 
nearly symmetric halves above and below the beam plane in the 
two regions of reduced charged-particle density which result 
from the sweeping action of the GOLIATH magnet. More de- 
tails about the experimental setup and the photon spectrometer 
can be found in 0] . 

The acceptance of the photon spectrometer for 7r° detec- 
tion in rapidity and transverse mass is shown in Fig. ^ The 
acceptance covers the region 2.3 < y < 3.0, near mid-rapidity 
(yen = 2.9). 

The minimum bias trigger requires a valid signal of the 
beam counters and a minimum amount of transverse energy 
5 GeV, detected by the Mid-Rapidity Calorimeter, 



> 



MIRAC piQ, which is located 24 m downstream of the tar- 
get. MIRAC consists of a hadronic and an electromagnetic sec- 
tion and covers the pseudorapidity interval 3.5 < ?/ < 5.5. 
Data have been taken with the 158 j4 GeV lead beam on tar- 
gets of Pb (495 and 239 mg/cm^) and Nb (218 mg/cm^). 
For the present analysis 9.7 million Pb-nPb and 0.23 million 
PbH-Nb minimum bias events were accumulated. The minimum 
bias cross sections have been calculated from the number of 
beam triggers and minimum bias triggers and the target thick- 
nesses. The yields have been corrected for small non-target 
background contributions (typically a few percent) to obtain 
(Tmb w 6300 mb and 4400 mb for Pb-nPb and PbH-Nb reac- 
tions, respectively. These absolute cross sections are estimated 
to have an overall systematic error of less than 10%. 



3 Pion Reconstruction and Efficiency 

The method for pion reconstruction is similar to that used by 
the WA80 collaboration as discussed in la]. It is also discussed 



M.M. Aggarwal et al.: Transverse mass distributions of neutral pions from Pb-induced reactions at 158-^ GeV 



3 



experimental setup of the 

WA98 experiment 
158 AGeV Pb+Pb collisions 
at the CERN SPS 



TOF-Wall (#2) 

PID of positive 
hadrons) 



Streamer-Tubes 
Pad-Chambers 



Goliath-Magnet 



Target 

(inside Plastic-Ball) 



Start- 
Counter 



highly segmented Lead Glass Calorimeter (LEDA) 

(identification of photons, Ji"- and rj-mesons) 



Forward 
Calorimeter 





Had.-Calorimeter 

(transverse energy) 



highly segmented Photon- 
Multiplicity Detector 



Veto-Detector 

for charged particles 



r TOF-Wall (#1) 

\^ (PID of negative hadrons) 



Multistep-Avalanche-Chambers 

with CCD-readout (track reconstruction 
of charged particles) 



. He in Silicon-Pad and 

the Target Region) Silicon-Drift Detectors 

(pseudorapidity distributions 
of charged particles) 



Fig. 1. The WA98 experimental setup. 



at length in [ ^5|^ ] so we will only briefly sketch the basics of 
the method. 

The showers used for the extraction of the neutral pion yield 
via the 77 decay branch can be selected with different criteria. 
It has been demonstrated that the extracted yield does not de- 
pend significantly on the choice of these photon identification 
cuts [P5|]. For the present analysis, showers with a small lateral 
dispersion have been selected as photon candidates. 

Hits in LEDA are combined pairwise to provide distribu- 
tions of pair mass vs. ttit — mo (where uit — \/pt + 
is the transverse mass and mp is the vr" rest mass) for all 
possible combinations. These distributions are obtained both 
for real events, R{minv, Jtit), and for so-called mixed events, 
M{minv, rriT), where a hit from one event is combined with a 
hit from another event with similar multiplicity. A/(TOinv, mx) 
provides a good description of the combinatorial background 
and is subtracted from R{minv, mr) to obtain the mass distri- 
bution of neutral pions (see Fig. ||). 

The large multiplicities, especially in central reactions of 
Pb -I- Pb, lead to a considerable probability that showers in the 
detector overlap and influence each other. Particles may be lost 
for reconstruction. In other cases the particle might be mea- 
sured with an incorrect energy. These effects lead to a detec- 
tor efficiency for pion reconstruction which depends on parti- 
cle density. To study this detection efficiency, the GEANT 
simulation package has been used to create artificial signals for 



the lead glass modules corresponding to neutral pions incident 
on the detector. The effects of detector noise and the digitiza- 
tion of the photomultiplier signals are also implemented. To 
obtain the 7r° reconstruction efficiency, these simulated pho- 
ton shower pairs have been superimposed onto the measured 
events. The simulation provides the means to extract the prob- 
ability that a pion at a given input transverse mass m^"* will be 
measured with urSp . These probabilities, extracted as a func- 
tion of 77i^\ can be understood as a matrix which transforms 
the real physical distributions into the measured ones. These 
correction matrices are used to correct the measured distribu- 
tions. A detailed description of the efficiency correction can be 
found in 

The distributions have then been corrected for contributions 
from reactions of Pb projectiles with material other than the 
target (e.g. residual gas, exit windows, etc.). The corresponding 
corrections were obtained by measurements performed without 
target. The corrections are negligible for medium-central and 
central reactions. 

The systematic errors on the measured transverse mass spec- 
tra are dominated by the following contributions: 

- An uncertainty in the absolute calibration of the momen- 
tum scale of 1 %. This may be translated into an uncertainty 
in the yields of a few % at low my rising to k, 13% at 
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Fig. 2. Top: Geometrical acceptance for neutral pions as a function of 
Pt and y. Bottom: Rapidity-integrated acceptance in 2.0 < y < 3.2 
as a function of tut — mo. 



- Uncertainties in the correction for contributions from no- 
target reactions, which are only relevant for peripheral re- 
actions. This leads to an error below 10% for ttit — toq > 
0.5 GeV/c2. 

The total centrality dependent systematic error is below 10% 
for rriT — rriQ > 0.5GeV/c^ for most centrality classes and 
below 20% for the most central sample. 



4 Results 

4.1 Neutral Pion Spectra 

Neutral pion spectra for Pb + Pb and Pb + Nb collisions for 
minimum bias and selected centralities are presented in Fig. ^ 
On the left hand side results for Pb + Pb minimum bias as well 
as for the 12.7% most central and the 17% most peripheral 
event selections are shown. Included also are results of expo- 
nential fits: 



/(rriT - mo) = C • exp 



rriT — mo 
T 



(1) 



lll-p-lllQ C 

0.9 - l.OGeV/c- 




0.2 0.3 

(GeV/c2) 



Fig. 3. Invariant mass distributions of photon pairs for minimum bias 
Pb -I- Pb collisions in different intervals of rriT — mo. The combina- 
torial background has been calculated using mixed-event distributions 
and subtracted. 



rriT — Too = 3.5 GeV/c^, which is independent of central- 
ity or reaction system. 

Uncertainties in the tt" extraction which include the er- 
ror in the determination of the invariant mass peak con- 
tent and the propagation of the energy resolution through 
the acceptance and efficiency corrections. This is largest 
in central collisions. It leads to an error below 10% for 
rriT — mo > G.5GeV/c^ for most samples. 



fitted over the range 0.7 GeV/c^ < mr - mo < 1.9 GeV/c^. 
As noted previously [|6|p7|], the data clearly deviate from the 
exponential shape when considered over the full range of trans- 
verse masses. The fit in the limited transverse mass range can, 
however, still be used to extract slope parameters T. For the 
12.7% most central collisions one finds T = 221 ± 2MeV 
which is significantly larger than the value of T = 208±5 MeV 
for peripheral collisions. Spectra for reactions of Pb + Nb (min- 
imum bias and 10% central) are shown on the right hand side. 
The spectral shapes are very similar to those for Pb + Pb col- 
lisions. For central collisions a slope parameter of T = 213 ± 
13 MeV is extracted which is intermediate compared to that of 
peripheral and central Pb + Pb collisions. Since the statistics 
for the Pb + Nb dataset is limited, we will concentrate on the 
Pb + Pb reactions in the rest of the paper 

In Fig. H the multiplicity distribution for central collisions 
is compared to predictions of the event generators FRITIOF 
7.03 @, VENUS 4.12 and HIJING 1.36 Clearly 
FRITIOF does not describe the data at all, while VENUS and 
HIJING yield a more reasonable description. The prediction 
of the pQCD calculation from [|l8[ is included as a solid line, 
which also shows a reasonable agreement with the data. The 
degree of agreement of the models can be better seen in Fig. ^ 
where the ratio of the experimental data to the generator re- 
sults is shown. FRITIOF is not included as the discrepancies 
are already evident from Fig. ||. VENUS overpredicts the pion 
production at high my by about a factor of two. HIJING shows 
the best agreement at large momenta but underpredicts the data 
slightly at intermediate my. Also the Monte Carlo models lead 
to a stronger concave curvature than the data. The pQCD calcu- 
lation overpredicts the data by w 30% in the range 1 . 5 GeV/c < 
mp — mo < 3.0 GeV/c. 

The momentum spectra for central collisions have been pub- 
lished akeady in JIt]]. As discussed in the corresponding erra- 
tum, cross section estimates had turned out to be incorrect. In 
the course of the present analysis it was realized that for the 
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Fig. 4. Invariant cross sections of neutral pions for Pb + Pb (left) and Pb + Nb collisions (right) of different centralities as a function of mr — mo. 



presentation of the results in figure 1 of Ref. [|l7|], the 7r° multi- 
plicities were incorrectly normalized. All other results and con- 
clusions of Ref. Jl^ ] are unchanged. Unfortunately, the same 
improperly normalized distributions were analyzed in Ref. mm . 
The correct normalization reduces the 7r° mulitiplicities by 27% 
with the result that the pQCD calculations presented there would 
overpredict the measured WA98 7r° result, as shown in Fig. ^. 
Still different parameters for the py-broadening in these calcu- 
lations might possibly enable the model to describe the central 
distributions better. More stringent tests of such a model can be 
performed when looking at the detailed centrality dependence 
of the pion production which has been done in the following. 

The minimum bias sample has been subdivided into eight 
centrality samples summarized in Table |l[ Measurements were 
performed with and without magnetic field, which does not af- 
fect the neutral pion distributions, but alters the transverse en- 
ergy used to determine the centrality. The corresponding cuts 
have been adjusted so that always the same fractions of the 
minimum bias cross section were selected. The data samples 
with and without magnetic field agree well with each other and 
have been combined in the present analysis. Table |l] shows the 
transverse energy cuts for one particular data set as an example. 

All the transverse mass distributions for Pb + Pb collisions 
for the different centrality classes as a function of mx — jtiq 



have a very similar shape (see Fig. 0). The distributions are 
compared more closely in Fig. ^ where the spectra for given 
centrality classes are divided by the minimum bias spectrum. 
Still in these ratios no drastic variations are seen when dis- 
played in logarithmic scale. The spectral shapes, especially for 
the intermediate centrality classes, are very similar. The spectra 
are broadened when going from peripheral reactions to semi- 
peripheral reactions, consistent with the nuclear enhancement 
seen in p H- A [Q] and in S H- Au collisions [|[|. However, the 
enhancement does not continue to grow for central collisions. 
Instead there is an indication of a stronger fall off in the most 
central class. 

This observation is surprising, because both initial state (par- 
ton multiple scattering) and final state (hadron rescattering) 
mechanisms are expected to yield a further broadening of the 
spectra with decreasing impact parameter, i.e. increasing thick- 
ness or volume. 

The change in the spectrum for the most stringent central- 
ity selection of 1% of the minimum bias cross section appears 
to be relatively strong, and the question arises whether this cut 
is really significant enough and whether this sample should not 
behave similarly to the adjacent sample. We have therefore per- 
formed detailed investigations of possible backgrounds and bi- 
ases with respect to this particular sample. Pile-up of multi- 
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pie beam interactions is effectively suppressed by strong time 
and amplitude cuts on the trigger level. In addition, the total 
calorimeter coverage of MIRAC and the ZDC is sufficient to 
reject all possible pile-up events, as the measured energy would 
exceed the beam energy. 

Also, simulations of the centrality selections including re- 
altistic fluctuations in the detectors show that e.g. classes 7 and 
8 are significantly different. This can e.g. be expressed by the 
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Fig. 7. Invariant multiplicities of neutral pions for Pb + Pb collisions 
of different centralities as a function of mr — mo. The solid lines 
show exponential fits to the spectra. 



mean (rms) of the distribution of the number of participating 
nucleons which are 346 (20) and 380 (10), respectively. 

As the centrality selection is performed with the transverse 
energy measured slightly off midrapidity, one might suspect 
that midrapidity particle production could suffer a different bias, 
if e.g. the pseudorapidity distribution would change signifi- 
cantly over the small range relevant here. This is, however, not 
observed. The charged particle pseudorapidity density at ry = 
increases in a fashion identical to the transverse energy as can 
be seen in Ref . [|T]] . 

While this leads us to expect no additional bias from the 
event selection, the direct comparison of the spectra obtained 
with different analysis cuts for the most central sample (class 
8) shows slightly larger variations than the other samples. We 
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have therefore assigned a larger systematic error of 20% to this 
sample, as was discussed in Section]^. 



class 


Et (GeV) 


0"/0"mb 


{Npart) 




1 


< 24.35 


17.2 % 


12 ±2 


9.9 ±2.5 


2 


24.35 - 55.45 


15.8 % 


30 ±2 


30 ±5 


3 


55.45 -114.85 


18.2 % 


63 ±2 


78 ± 12 


4 


114.85 -237.35 


23.5 % 


132 ±3 


207 ± 21 


5 


237.35 -326.05 


12.3 % 


224 ± 1 


408 ± 41 


6 


326.05 -380.35 


6.2 % 


290 ±2 


569 ± 57 


7 


380.35 -443.20 


5.8 % 


346 ±1 


712 ±71 


8 


>443.20 


1.0 % 


380 ±1 


807 ±81 


6-8 


>326.05 


12.7 % 


323 ± 1 


651 ± 65 



Table 1. Centrality classes, as selected by the amount of transverse 
energy measured in MIRAC, for Pb -I- Pb collisions. The average num- 
ber of participants and binary nucleon-nucleon collisions as calculated 
with VENUS 4. 12 with an estimate of the systematic error. Please note 
that the systematic error in the number of collisions is correlated for 
all samples. 
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Fig. 9. Inverse slope parameters T of neutral pion transverse mass dis- 
tributions as a function of the number of participants in Pb -I- Pb colli- 
sions. The exponential fits are performed in different regions rriT — mo 
as indicated. The error bars shown contain statistical and systematic 
errors added in quadrature. 



4.2 Average and Inverse Slopes 

Another means to characterize the spectral shape is via the av- 
erage transverse momentum {pr), or via the inverse slopes 
T of the spectra. The truncated av erag e transverse momen- 
tum has already been presented in [|l7|], where it was shown 
that the values of (pt(Pt""))' ^^'-^ ^ lower cutoff of p™"' = 
0.4 GeV/c^, increase from peripheral to medium central colli- 
sions but seem to saturate for still smaller impact parameters. 
Since the (pt(p™™)) is dominated by the momentum region 
near the cutoff, this is not in contradiction with the dependence 



observed in Fig. [s] A similar analysis has been performed with 
inverse slope parameters which were obtained by fitting expo- 
nentials (Eq. I]) to the spectra in limited regions of transverse 
mass. The extracted slope parameters as a function of the num- 
ber of participants are shown in Fig. ^ Fits in the lowest rriT 
interval yield a slope of T = 204 McV for peripheral reac- 
tions. In the intermediate centrality range the slope appears to 
be constant at T w 220 MeV. 

Inverse slopes have also been extracted for other inter- 
vals of transverse mass (see Fig. ||). For all centralities the 
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Fig. 10. Exponent a of the mass dependence of neutral pion cross sec- 
tions as a function of transverse momentum. Values are obtained from 
ratios of minimum bias Pb + Pb and Pb + Nb collisions to a param- 
eterization of p-l-p collisions. Included for comparison are results for 
S-l-Au collisions at 200 AGeV and the proton induced data from 
Ref. [§]. 



slopes increase with larger transverse mass, which is another 
demonstration of the curvature of the spectra. The centrality 
dependence is more pronounced for the higher uit regions 
(1 — 2 GeV/c^ and 2 — 3 GeV/c^). There is a continuous rise in 
the inverse slope from very peripheral reactions up to reactions 
with Npart ~ 130, the highest slopes reaching T « 240 MeV 
and T w 260 MeV, respectively. There is an indication of a de- 
crease in the slope for very central collisions which is however 
not conclusive in view of the systematic errors. 

From this analysis of the spectral shapes by looking at spec- 
tral ratios or inverse slopes one may conclude, that apart from 
the broadening in going from peripheral to semi-peripheral 
collisions, no striking features are observed. Apparently these 
analysis tools are not suited to extract information on possible 
more subtle variations. 



4.3 Scaling with System Size 

In addition to variations of the shape of the momentum spectra 
it is of interest to study the variations in absolute multiplicities. 
Especially at high transverse momentum one naively expects 
an increase of the cross section proportional to the mass num- 
ber of the nuclei and, correspondingly, an increase of the mul- 
tiplicity proportional to the number of binary collisions due to 
the importance of hard scattering. In fact, it was already ob- 
served in pA collisions at beam energies of 200 — 400 GeV [Q] 
that the increase in cross section at high transverse momenta is 
even stronger than the increase in the target mass. 

This behaviour can be investigated by parameterizing the 
data according to the phenomenological expression: 



PJ 



(2) 



where A and B are the mass numbers of the projectile and tar- 
get nuclei. For this purpose the experimental data have been 
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Fig. 11. Exponent a of neutral pion transverse mass distributions pa- 
rameterized as a function of the number of participants (a) and colli- 
sions (b) in Pb + Pb collisions. The exponents are shown as a function 
of pt- The open circles are obtained from the ratio of semi-peripheral 
Pb+Pb data to p+p data, while for the filled circles fits of a power law 
to the seven most central classes are performed. The grey bands indi- 
cate the systematic error relevant for this fit from the calculation of the 
number of participants and collisions and from the centrality depen- 
dent systematic error of the neutral pion multiplicity. The error bars 
are only statistical. 



compared to a parameterization [ [26[ ] of pp reactions [ |32[p3 34[p5y3q , |37| 
scaled to the same y^: 



E- 



C ■ 



Po 



(3) 



dp^ \Pt+Po 

with C = 4.125c3/GeV^ pa = 9.02 GeV/c and n ^ 55.77. 
We have assigned a systematic error of 20% to this reference 
distribution. The values of the exponents a for minimum bias 
reactions of Pb+Pb and Pb+Nb are shown in figure |l^ Included 
for comparison are results for Sh-Au collisions at 200 j4GeV 
and the proton induced data from [Q]. All data follow the 
same trend, i.e. the exponent is considerably lower than one 
for low transverse momenta and increases monotonically with 
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normalized to the number of binary collisions, also called the nuclear 
modification factor. The grey band shows the estimate of the system- 
atic error due to the calculation of the number of collisions and the 
absolute cross section normalization relative to p+p. 



Pt finally reaching values > 1 for highpT- However, the val- 
ues differ by about 10% for the different systems and the point 
where they cross the line a = 1 ranges frompT = 1-5 GeV/c^ 
to 2.5 GeV/c^. It should be noted that the different experimen- 
tal trigger biases may influence these results. 

This form of scaling analysis is only applicable to mini- 
mum bias reactions. To study the scaling behaviour for heavy 
ion reactions of different centralities the scaling with the num- 
ber of participants or nucleon-nucleon collisions is more ap- 
propriate. An analysis of the scaling of the charged multiplic- 
ity with the number of participants has been performed in [Sop . 
A more general analysis of the scaling behaviour of transverse 
energy and particle production in Pb + Pb collisions has been 
presented in ^)\\ and the scaling of neutral pion production 
with the number of participants in [p^. For the present paper 
the number of participants and collisions have been obtained as 
described in [plt|. The values are given in Table |l[ The system- 
atic errors of these quantities have been obtained by comparing 
calculations with different assumptions. In addition to VENUS 
and FRITIOF calculations with default parameter settings, we 
have investigated VENUS calculations with a modified nuclear 
density distribution, with a significantly worsened energy reso- 
lution of the calorimeter, with different assumptions on the size 
of the minimum bias cross section and ignoring the target-out 
contribution. These calculations are described in pi]]. More- 
over we have also investigated independently results of a toy 
model calculation of nuclear geometry and also the influence 
of modest variations of the nucleon-nucleon cross section en- 
tering into the calculations. The maximum deviation of any of 



these other calculations from the default one is within the error 
given in Table |l]. As expected, the relative error is largest in the 
most peripheral sample. This is due to the uncertainty in mod- 
elling the i?T-fluctuations and the trigger threshold — effects 
which have been included in our error estimate. 

It should be noted that the most peripheral data set in this 
analysis does not correspond to the most peripheral collisions 
theoretically possible from nuclear geometry, as they are re- 
jected by the minimum bias trigger. Our minimum bias trig- 
ger cross section of « 6300 mb corresponds to w 85% of the 
expected total geometrical cross section. The more peripheral 
reactions rejected by our trigger would of course have a still 
larger uncertainty due to e.g. the uncertainties in the nuclear 
density distribution. In our samples the error is still at a reason- 
able level as quoted in the table. 

Furthermore, the error in ratios of e.g. the number of colli- 
sions for different samples is considerably smaller, since some 
of the systematic effects are correlated, as is obvious e.g. for 
the size of the nucleon-nucleon cross section. 

The scaling of the neutral pion production is analyzed both 
as a function of the number of collisions and as function of the 
number of participants. As in the earlier publications the yield 
is parameterized as: 



(Npart) OC (Npart) 



or 



(4) 



(5) 



The results are presented in Fig. [IT] where the upper panel 
shows the scaling with the number of participants and the lower 
panel with the number of collisions. The open circles represent 
the exponents calculated from the ratio of semi-peripheral col- 
lisions (class 2) to the parameterization of pp collisions (Eq. |^). 
In this case a behaviour very similar to the original Cronin- 
effect is observed: The exponents increase monotonically with 
increasing pT reaching values > 1 for large px- A scaling of the 
particle yields with (Ncou)""^^ is equivalent to a scaling of the 
minimum bias cross sections with A"^^. This value is reached 
in Fig. [rT| b at px ~ 1.5 GcV/c^, similar to where the scaling 
exponents of the p+A and S+Au data in Fig. |l^ reach the value 
of 1 . However, the minimum bias Pb+Pb data in Fig. |o| show 
a weaker scaling exponent at intermediate and higher px- 

Fits to the centrality classes 2-8 (displayed as filled circles 
in Fig. [ri|), i.e. up to the most central collisions, show almost 
constant exponents for all transverse momenta. In fact, the val- 
ues even indicate a slight decrease with increasing px- At the 
highest px the exponents appear to be ctp w 1 and ac w 0.8. 
They are considerably below the scaling exponents for semi- 
peripheral Pb+Pb relative to p+p. Noticeably the scaling expo- 
nent ac is also significantly below the value of 1 expected as a 
naive scaling for hard scatterings. 

Ratios of the measured pion multiplicity distributions for 
two different samples (labeled X and Y) normalized to the num- 
ber of collisions given in Table [l|: 



Rxvimx) 



E 



'^{mx)/NcoU 



X 



E 



{mx)/Ncoii 



(6) 
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Fig. 13. Ratios of invariant multiplicity distributions of neutral pions 
for central to peripheral Pb+Pb collisions normalized to the number 
of binary collisions. The left plot shows the ratio using the most pe- 
ripheral sample, the right plot a similar ratio with the second most 
peripheral sample. The grey bands show the estimate of the system- 
atic error due to the calculation of the number of collisions and due to 
the systematic error in the tt" distribution. 



are shown in Figs. |l2|, |l3] and [ij. In Fig. [l^ the ratio of the 
12.7% most central collisions to the parameterization of p+p is 
shown: 



E 



Rximr) 



{mT)/Ncoii 



X 



Td^N 



{mr) 



(7) 



This special case of equatio n N is sometimes referred to as the 
nuclear modification factor [|4^. The ratio is k, 0.3 at low rriT 
and increases exponentially towards higher rriT approaching 
values close to 10. [] This is in line with the earlier observations 
and demonstrates again that relative to p+p reactions there is 
a strong Cronin effect in heavy ion collisions for all centrali- 
ties. However, as evident from Fig. |ll]b, most of the anomalous 
enhancement is already present in peripheral reactions, and the 
additional enhanement in central reactions is much weaker. In 
the following we will concentrate on this particular observa- 
tion. 

Fig. |lj shows ratios of the 12.7% most central collisions 
to the two most peripheral samples. Both ratios are below 1 
for low niT- In comparison to the most peripheral sample (left 



' The same ratio for central Pb+Pb(Au) collisions is shown in |^]. 
There a curve is drawn which saturates at i? = 2 for high px- It is ac- 
tually dominated by data points at pt = 2 — 3 GeV / c, and the WA98 
data points can be seen to lie above the curve at high pr- Remaining 
differences in the two figures might be due to different estimates for 
the pp-spectra being used in the ratios. 
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Fig. 14. Ratios of invariant multiplicity distributions of neutral pions 
normalized to the number of binary collisions. The upper left plot (a) 
shows the ratio of peripheral Pb+Pb collisions to the parameterization 
of p+p [^]. The other plots (b-d) show different ratios of a more central 
sample to a more peripheral sample of Pb+Pb collisions. The filled cir- 
cles show the experimental results, the open circles are results from the 
HlJlNG event generator and the solid lines similar ratios from pQCD 
calculations including pT broadening The grey bands show the 
estimate of the systematic error due to the systematic error in the tt" 
distribution. The additional error bar in (a) and (b) shows the normal- 
ization uncertainty due to the calculation of the number of collisions 
and the absolute cross section normalization relative to p+p. 



part) the spectra are compatible with a scaling with the num- 
ber of collisions in central reactions at higher uit, while rela- 
tive to the semi-peripheral class (right part) the enhancement is 
significantly weaker even at high m^. Similar ratios for more 
detailed centrality selections are shown in Fig. |l4|. T he ratio of 
peripheral Pb+Pb collisions to p+p (Eq. H) (Fig. |14^) increases 
strongly with increasing transverse mass - this is in line with 
the Cronin effect discussed above. A similar trend is observed 



when going from peripheral to medium-central data (Fig. 143). 
In addition, the pion production is seen to increase roughly pro- 
portional to the number of collisions even at low transverse 
masses. Going from medium central to central (Fig. |l^) the 
trend is reversed: the ratio decreases with increasing transverse 
mass as was already seen in Fig. ^ and the pion multiplicities 
increase more weakly than the number of collisions. The ratio 
of very central to central collisions shows an indication of a 
similar effect although not very significant. 

Included in Fig. [l^are results of HIJING calculations (open 
circles) for the same centralities. They show a very different 
trend: For all but the most central case the ratios increase with 
increasing transverse mass, for high transverse masses the ra- 
tio is always larger than one. HIJING does not describe p+p 
data well at these energies while a reasonable description of 
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Fig. 15. Multiplicities of neutral pions normalized to the number of 
binary collisions Run = {Ed^ N / dp^) / Ncoii as a function of Ncoii 
for different ranges of transverse mass as indicated in the pictures. The 
filled circles show the experimental results for Pb+Pb reactions, the 
open circles the parameterization for p+p reactions (Eq. ^. The error 
bars contain statistical and systematic errors added in quadrature. 



central Pb+Pb collisions is obtained by an implementation of 
the Cronin effect via a soft transverse momentum kick model 
which introduces a very strong A dependence, as was pointed 
out already in [^|. Thus the reasonable description of central 
Pb+Pb collisions by HIJING appears to be fortuitous. The au- 



thors of 1 41 1 stated that a better implementation of the Cronin 
effect should use an additional intrinsic pT of the incoming par- 
tons and a likely weaker pT broadening. Still, in any model of 
the Cronin effect one would expect the ratios of central to less 
central spectra, as shown in Fig. p^j s, to be larger than one at 
high transverse mass. Similarly, one would expect the scaling 
exponent shown in Fig. Illjto be ac > 1 at high transverse mo- 
mentum. A more refined calculation including intrinsic pT was 
performed in using the same model as in JT^. The results 
are also shown in Fig. |l^ as solid lines. As expected, the ratios 
show a weaker increase at high riiT compared to HIJING, but 
are all > 1 and thus do not explain the centrality dependence 
as seen in Fig. |l4| :. 

With these observations in mind one can revisit the evolu- 
tion of the particle yields with system size at given transverse 
mass. This is done in Fig. [l^ where the neutral pion yield per 
binary collision; 



bin 



(8) 



is shown as a function of the number of collisions for four dif- 
ferent transverse mass intervals. In addition to Pb+Pb reactions 
also the parameterization of p+p is included. At relatively low 
transverse mass the yield per collision decreases as expected 
from the scaling exponent ac < I- With increasing transverse 



mass a rise of Rbm develops which is most prominent for the 
highest rriT- However, it is also clear that i?;,i„ decreases signif- 
icantly with Ncou for Ncou > 200. The exact point of turnover 
of the evolution of the pion yields could not be so easily ob- 
served in the investigation of the scaling exponents above - it 
is much more clearly born out when normalizing to the number 
of binary collisions as in Rbm- 



5 Summary 

We have presented transverse mass spectra of neutral pions 
in Pb+Pb reactions at 158-^ GeV for different centralities. A 
comparison to several event generators found none of them able 
to adequately describe the spectra. 

While there is a strong broadening of the transverse mass 
spectra in going from p+p to peripheral Pb+Pb reactions, and 
further to semi-peripheral reactions, there is an indication of 
a stronger fall off of the spectra when going to the most cen- 
tral selections. The absolute yields at high transverse momenta 
show an enhancement which grows stronger than with the num- 
ber of collisions up to medium central Pb+Pb collisions, which 
is qualitatively similar to the Cronin effect observed in p+A 
collisions. For central collisions, however, the further increase 
in multiplicity at high pT is weaker than with the number of 
collisions. This is in qualitative contradiction to conventional 
explanations of the Cronin effect (py-broadening) which is ex- 
pected to cause a further strengthening of the nuclear enhance- 
ment with more central collisions. As a possible explanation 
of this behaviour one may consider that the multiple scatter- 
ing mechanisms which are expected to be responsible for the 
apparent p^-broadening in pA and peripheral Pb+Pb might be 
modified in central Pb+Pb collisions. This might be possible 
if both initial and final state scatterings are relevant for the 
enhancement, and the relative contributions are shifted more 
towards final state contributions in central collisions, reminis- 
cent of a more and more thermalized system. Alternative expla- 
nations might involve suppression mechanisms independent of 
the nuclear enhancement in question, e.g. an onset of quench- 
ing via energy loss of produced particles (partons or hadrons) 
in central Pb+Pb collisions. 
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